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Xn  recant  years,  tree-structured  analy- 
ala/reconstruction  ayatema  have  been  extensively 
studied  for  use  in  subband  coders  for  speech,  la 
such  systems,  it  is  important  that  the  individual 
channel  signals  be  decimated  in  such  a  way  that 
the  number  of  samples  coded  and  transmitted  does 
not  exceed  the  number  of  samples  in  the  original 
speech  signal.  Under  this  constraint,  the 
systems  presented  in  the  past  have  sought  to 
remove  the  aliasing  distortion  while  minimising 
the  overall  analysis/reconstruction  distortion. 
Xn  this  paper,  it  is  shown  that  it  is  possible  to 
design  tree-structured  analysis/reconstruction 
systems  which  meet  the  sampling  rate  condition 
and  which  alao  result  in  exset  reconstruction  of 
the  input  signal.  This  paper  develops  the  condi¬ 
tions  for  exact  reconstruction  and  presents  a 
general  method  for  designing  the  corresponding 
high  quality  analysis  and  reconstruction  filters. 


Xn  a  recent  paper  by  Barnwell  [1],  it  was 
shown  that  alias-free  reconstruction  using  recur¬ 
sive  and  non-recuralve  filters  was  possible  where 
the  analyilt/reconstructlon  section  had  no  fre¬ 
quency  distortion  or  no  phase  distortion,  but  not 
both.  Xn  the  development  that  follows,  the  coef¬ 
ficient  spmttry  condition  on  the  analysis 
filters  is  lifted  end  exact  reconstruction  free 
of  aliasing,  phase  distortion  and  frequency  dis¬ 
tortion  is  shown  to  be  possible  using  FIR 
filters.  Xn  addition,  the  filter  constraints 
that  enable  perfect  reconstruction  are  discussed 
and  an  easily  implement able  design  procedure 
providing  high  quality  filters  is  presented. 


The  frequency  division  of  the  subband  coder 
is  performed  in  the  analysis  stgoe  as  shown  in 
Figure  1,*  where  B  (e3w)  end  B.U7*)  are  lowpaas 
and  hlghpasa  filters  respectively.  To  preserve 
the  system  sapling  rate,  both  channels  are 
decimated  resulting  in  the  two  down-aanpled 
signals,  Y  <e'*)  and  Y.<eO.  In  the  recon¬ 
struction  action,  the  Sands  are  recombined  by 
up-sopling  end  filtering  to  give  the  recon¬ 
structed  signal 


£<«3“>“j[ho(.3*i)go(.3“>«s1(.3*i)g1<.3,“)  )«(.3“) 

(D 

The  frequency  response  of  the  2-band  system  is 
represented  by  the  first  term  in  equation  (I), 
while  the  second  tern  is  the  aliasing* 

TO  obtain  exact  reconstruction,  consider  the 
esse  in  which  the  analysi e/reconstruct Ion  filters 
ere  designed  using  the  assignments 

C(.3“)-I(.”1*‘)  (3.) 

o  o 

C,  (2b) 

H,  <«3tl)-B0<-."3“)  (2c) 

where  B  (e  is  not  required  to  be  linear  phase, 
o 

This  assignment  eliminates  the  aliasing  term 
in  equation  (1)  and  the  resulting  overall  system 
function,  is  given  by 

C<ej  V 

•  (3) 

where  F  (e^  end  F.fe^  will  be  called  the 
"product  filters."  *Xt  la  clear  from  equation  (3) 
that 


*1<«3'“>-ro<-.34') 


Bow  assume  that  the  product  filters  are  sero 
phase  FXR  filters.  Then  the  exact  reconstruction 
condition  is  given  by 

c(n)  *  {C(e^  }  ■  4(a)  (4) 

where  denotes  the  inverse  Fourier  trans¬ 
form.  The  corresponding  condition  on  the  product 
filter  is  glvsn  by 


Vn>  *  «(•> 
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So  In  order  for  exact  reconstruction  to  be 
obtained ,  there  arc  two  conditions  on  the  product 
filters  which  oust  be  met.  first,  they  oust  seet 
the  condition  of  equation  (5).  Second,  they  oust 
be  decomposable  Into  analysis  and  reconstruction 
filters  Ip  such  a  way  that  equation  (3)  is  valid. 

Xn  order  to  see  how  these  conditions  can  be 
Mt,  It  is  appropriate  to  decompose  the  product 
filters  into  their  sero  time  and  nonsero  time 
components,  giving 

(«) 

in  the  frequency  domain  or,  equivalently, 

f0(n)-v(n)4A5(n)  (7) 

In  the  time  domain.  Xn  these  expressions,  v(n) 
is  constrained  to  have  no  sero  time  component, 
i.e.»  v(0)  •  0,  and  the  exact  reconstruction 

condition  of  equation  (10)  is  met  so  long  as  A  * 
1  and 

v(2n)4)  n  •  0,*,1,*,2, ...  (•) 

or,  equivalently 

<*} 

This  la  the  class  of  all  filters  which  are  anti¬ 
symmetric  about  the  Nyqulst  frequency,  as 
illustrated  in  Figure  3.  This  constraint  can  be 
met  In  several  ways  and  a  number  of  powerful 
filter  design  tools  are  available  to  design  tbe 
required  filters. 

Filter  Peal— 

Xn  the  design  procedure,  the  approeeh  is  to 
first  design  the  product  filter,  F^*)*  end  then 
to  decompose  F0(s)  into  B^U)  end  BpU*1).  there 
are  e  number  of  techniques  which  can  be  used  to 
design  Fe(s). 

For  eaumpls,  any  filter  of  the  form 
sin  r 

f  (n)»w(n) - — >  (10) 

SB 


where  w(n)  is  a  window  function,  will  satisfy  the 
exact  reconstruction  condition  of  equation  (S). 
Similarly,  optimal  equi ripple  product  filters  can 
be  constructed  using  a  Asset  exchange 
algorithm.  Xn  such  designs  (in  which  the 
Chebysbev  error  is  minimised  with  equal  weighting 
over  the  entire  frequency  bend) ,  only  those  solu¬ 
tions  with  the  largest  number  of  extrema,  i.e. 
the  extraripple  filters,  will  satisfy  equations 
(•)  and  (t).  Both  the  Parka-McClellan  end  the 
■ofstetter  algorithms  for  equlrlpple  filter  de¬ 
sign  may  be  used  to  design  the  product  filters. 

The  attenuation  of  tbe  analysis  filter  and 
product  filter  ere  related  in  the  following  may 

af  *  ♦  aoiog10<a-io  •  )  (11) 


where  Ip  la  the  product  filter  attenuation  and  A. 
is  the  'analysis  filter  attenuation.  The  secono 
term  in  this  expression  is  dus  to  the  difference 
in  the  transition  width  of  FQ(s)  and  BQ(s).  Thia 
tends  to  make  the  Hofstetter  algorithm  (2)  more 
attractive  for  the  equal  ripple  designs  slnoa  tbs 
attenuation  la  an  explicit  input  to  the  design 
procedure.  Xt  is  important  to  notice  that  tbe 
transition  width  of  the  analysis  filter  will  be 
larger  than  that  of  the  product  filter  end  there¬ 
fore  care  should  be  exercised  in  using  the  Farts- 
McClellan  algorithm. 

The  second  condition  which  must  be  met  Is 
that  tbe  product  filters  must  be  factorable  in 
such  a  way  that  equation  (2c)  is  valid.  Tbe 
additional  constraint  that  equation  (2c)  places 
on  Fc(t)  can  be  stated  *s  follower  for  every 
sero  of  F0(t)  et  re3  ,  there  must  be  another 
corresponding  taro  at  (1/r)e3  .  Xf  thia  condi¬ 
tion  is  met,  then  FQ(z)  can  always  be  written  as 
tbe  product  of  reciprocally  paired  seroe 

*  -1 

F  (*)<  8  (s— *  Ms  A-S  )  (12) 

°  _  ■  • 


where  G  la  a  real  constant.  An  important  point 
to  note  la  that  any  FDR  filter  which  is  both  real 
and  whose  coefficients  are  symmetric  in  the  time 
domain  (sero  phase)  comes  close  to  meeting  this 
condition.  Xn  particular,  the  zeros  of  any  such 
real,  syne  trie  FXJt  filter  must  either  meet  the 
above  condition,  lie  in  complex  sero  pairs  on  the 
unit  circle  or  occur  et  a  •  el.  These  two  condi¬ 
tions  become  identical  if  there  Is  a  further 
requirement  that  any  sero  on  the  unit  circle  must 
he  a  double  sero.  The  Important  point  bars  is 
that  any  product  filter  which  meets  the  exact 
reconstruction  oondltlon  of  equation  (10)  can  he 
easily  transformed  into  e  new  product  filter 
which  also  meets  tbs  factor  isatlon  oondltlon 
required  to  satisfy  equation  (2c) •  Tbe  required 
transformation  la  simply 

f  (#)  ♦  af  <n)H>  (13) 

o  o 

where  "a*  and  "b"  are  real  constants. 

This  is  exactly  the  strategy  outlined  by 
Bsrrnann  end  tch ussier  (3)  for  designing  minima 
phase  FXA  filters.  The  starting  point  is  e 
symetcic,  sero  phase  half-band  filter  with  an 
odd  number  of  coefficients.  Figure  3  shows  e 
real  and  sysmetric  product  filter  (equation  13) 
with  *b*  equal  to  sero.  Picture  now,  if  you 
will,  the  earns  filter  with  "b"  equal  to  a  half. 
The  stopband  ripples  now  occur  above  and  below 
the  real  axle.  These  sero  crossings  mark  the 
location  of  roots  on  the  unit  circle.  As  "b"  is 
increased,  the  sywetry  conditions  ere  mot 
changed,  but  the  sere#  on  the  tat  it  circle  migrate 
toward  one  another  end  leave  tbe  unit  circle  in 
pairs.  Bben  "b*  is  large  enough,  all  of  the 
seroe  will  have  either  been  driven  off  the  unit 
circle  or  will  have  become  double  seroe  on  the 
unit  circle  (Figure  2).  For  equlrlpple  filters, 
it  is  possible  to  have  many  double  teroe  on  the 


unit  circle,  while  for  window  design#  there  will 
generally  only  be  one  double  zero  pair.  The 
effect  of  the  multiplier  constant  "a"  is  to  re- 
scale  the  Impulse  response  to  have  the  correct 
passband  gain.  It  can  be  easily  shown  that 


-A/20 

*.1+10  F 


(14a) 


-*  /20 

to-  p  (14b) 


Once  the  appropriate  product  filters  have 
been  designed,  it  is  straightforward  to  extract 
the  analysis  and  reconstruction  filters.  Per 
this  case,  the  product  filter  can  be  characteris¬ 
ed  as  in  equation  (12)  and,  therefore 

*  -1 

F.  (*)  -  F  (— x)  -  G  H  («*z  )(z  ♦*  )  (15) 

i  °  ml  ■  ■ 

where  G  is  a  constant.  The  corresponding  analy¬ 
sis  and  reconstruction  filters  are  given  by 

M 

Be(*)-/5“j|l(«-»>)  (ito) 


remaining  roots  comprise  the  corresponding  recon¬ 
struction  filter  which,  in  this  case,  would  have 
maximum  phase.  Figs.  4  and  5  show  the  normalised 
magnitude  response  and  group  delay  respectively 
of  a  32  coefficient  minimum  phase  analysis 
filter.  Also  shown  (Fig.  €)  is  an  approximate 
linear  phase  analysis  filter  obtained  by  select¬ 
ing  conjugate  pairs  of  product  filter  roots  that 
occur  alternately  inside  and  outside  of  the  salt 
circle  as  the  frequency  increases. 

Third,  note  that  the  analysis  filters  are 
related  in  that  h1  (n)  is  formed  by  reversing 
hpta)  in  time  and  multiplying  by  (-1)*.  the 
implication  is  that  the  filter  coefficient  symme¬ 
try  present  in  QKFs  is  absent  here  and  therefore 
the  polyphase  structure  (1,4)  cannot  be  used  in 
implementation. 

The  major  point  is  that  exact  reconstruction 
analysis  and  reconstruction  filters  are  not  dif¬ 
ficult  to  design.  In  addition,  the  exact 

reconstruction  filters  generated  by  this  techni¬ 
que  generally  have  better  characteristics  in 
terms  of  their  transition  widths  and  attenuation 
than  the  most  published  QKFs  of  the  same  length, 
and,  of  course,  they  also  give  exact  recon¬ 
struct  ion. 


s  u"l4*  > 
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Filters  of  this  type  give  the  desired  exact  re¬ 
construction. 
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A  close  examination  of  equations  (14a-1dd) 
reveals  some  interesting  points.  Firat^  each 
sero  pair  in  Fp(s)  (one  at  re3*  and  one 
at  (1/r)e3*)  one  of  the  two  zeros  is  always  in¬ 
cluded  in  the  analysis  filter  while  the  other  is 
Included  in  the  reconstruction  filter.  Hence, 
the  analysis  filters  and  the  reconstruction 
filters  always  have  identical  magnitude 
responses,  and  this  magnitude  response  is  exactly 
the  square  root  of  F  (e3w).  Therefore,  F  (e3*5 
must  be  designed  to  Be  the  square  of  the  desired 
analysis  filter  magnitude. 

Haoond,  note  that  for  a  given  product 
filter,  there  are  many  choices  of  >0(t)  that  have 
a  magnitude  response  equal  to  the  square  root  of 
the  product  filter  magnitude.  Among  these  is  the 
minimal  phsse  analysis  filter  consisting  of  roots 
strictly  Inside  the  unit  circle  end  one  root  from 
each  double  sero  pair  on  the  unit  circle.  Tbs 


(4]  R.  X.  Crochlere,  "On  the  Design  of  Bub-Band 
Coders  for  Low-Bit-Rate  Speech  Cornwall- 
cations,"  Bell  System  Tele.  Journal,  94, 
(May-June,  1977),  pp.  747-770. 

15]  A.  Croleier,  D.  Xsteban  and  C.  Gal  and, 
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FIGURE  1.  2-BAND  SUBBAND  CODER. 


FIGURE  2.  POLE/ZERO  PLOT  OF  PRODUCT 
FILTER 


FIGURES.  NONZERO  TIME  COMPONENT  OF  PRODUCT 
FILTER.  VI#  PI. 
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FIGURE  E  32-TAP  MINIMUM  PHASE  FILTER. 


FREQUENCY 

FIGURE  4.  32-TAP  MINIMUM  PHASE  FILTER. 


FIGURE  «.  32-TAP  APPROXIMATELY  LINEAR 
PHASE  FILTER. 


